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FOREWORD 


This  report  covers  the  progress  made  In  the  third  and  final  year, 

22  January  1973  through  2 February  1974,  on  a research  program  con- 
ducted by  the  Department  of  Aerospace  Engineering,  The  University  of 
Michigan,  Ann  Arbor,  Michigan,  under  Contract  F08635-71-C-0083 
with  the  Air  Force  Armament  Laboratory,  Armament  Development  Test  Center, 
Eglin  Air  Force  Base,  Florida.  Dr.  J.  C.  Foster,  Jr.  (DLJl)  served 
as  program  manager  for  the  Armament  Laboratory.  Contractor  research 
was  directed  by  Professor  J.  A.  Nlcholls. 

This  technical  report  has  been  reviewed  and  is  approved  for 
publication. 

FOR  TOE  COMANDER 


ABSTRACT 


A simplified  version  of  a combustible  cloud  is  considered  wherein  the 
cloud  is  homogeneous  throughout.  Blast  wave  initiation  of  detonation  in 
this  cloud  is  treated.  The  self  •'Similar  character  of  strong  blast  waves 
and  Chapman- Jouguet  detonation  waves  is  used  to  arrive  at  simplified 
closed  form  solutions  for  the  generation  of  ground  impulse  and  dynamic 
Impulse,  up  to  the  time  that  the  detonation  arrives  at  the  edge  of  the 
cloud.  As  a result, assessment  of  the  effects  of  fuel  properties  and  cloud 
geometry  are  readily  obtained.  The  influence  of  side  relief  at  the  top  of 
the  cloud  on  the  ground  impulse  is  also  considered. 

Current  experimental  research  is  described  which  investigated  the 
initiation,  transition  and  quasi-steady  propagation  processes  associated 
with  blast  Initiated,  cylindrical  detonation  waves.  Experiments  were 
conducted  with  gaseous  and  heterogeneous  fuel-air  mixtures,  using  a 
specially  designed  sectored  shock  tube,  and  the  results  compared  with 
simplified  analytical  descriptions.  Gas  phase  detonation  research  has 
provided  insight  into  both  overall  and  detailed  wave  behavior  for  a full 
volumetric  range  of  Methyl  Acetylene,  Propane,  Propadlene  (MAPP) -air 
mixtures.  Heterogeneous  detonation  research  has  allowed  examination 
of  gross  wave  behavior  for  limited  test  conditions  and  additionally  provided 
details  as  to  the  time  history  of  drop  breakup  and  ignition  when  impacted 


by  strong  blast  waves.  Both  of  these  phases  provide  evidence  to  substanti- 
ate, in  accordance  with  existing  approximate  theories,  that  the  strength 
of  the  resultant  combustion  wave  generated  is  dependent  upon  the  blast 
source  energy. 

Variations  in  fuel  concentration  and  blast  source  energy  were  found 
to  have  a predictable  effect  upon  wave  pressure,  propagation  velocity, 
detonation  limits,  Ignition  threshold  energy  limits,  wave  transition  radius, 
and  wave  transition  structure.  Comparisons  made  with  the  work  of  others 
reveal  these  findings  do  not  suffer  from  scale  effects. 
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SECTION  I 
INTRODUCTION 


Many  facets  of  unconfined  explosions,  such  as  a fuel-air 
explosion,  are  ill  understood.  Theoretical  treatments  as  well  as 
carefully  controlled  experimental  studies  of  the  problem  have  been 
relatively  few.  The  aim  of  this  research  program  is  to  add  to  the 
understanding  through  an  integrated  analytical  and  experimental 
study.  The  idea  has  been  to  treat  a simplified  model  of  the  ex- 
plosive cloud  so  that  general  conclusions  can  be  drawn.  For 
example,  it  is  desirable  to  understand  the  pressure  and  velocity 
fields  that  are  set  up,  the  impulse  characteristics,  the  influence  of 
fuel  properties,  the  influence  of  cloud  geometry  and  non-uniform 
distribution  of  fuel  throughout  the  cloud,  the  detonation  initiation 
requirements,  etc.  Many  of  these  subjects,  but  net  all,  are  treated 
in  this  report.  The  theoretical  aspects  are  taken  up  in  Section  II 
and  the  results  of  experimental  studies  are  discussed  in  Section  in. 


SECTION  n 

THE  IMPULSE  GENERATED  BY  BLAST  WAVES  PROPAGATING 
THROUGH  COMBUSTIBLE  MIXTURES 

A.  INTRODUCTION 

Detonations  can  be  initiated  in  combustible  mixtures  by  suf- 
ficiently strong  blast  waves.  Such  blast  initiated  detonations,  or  fuel- 
air  explosions,  are  of  interest  in  certain  military  applications,  and 
may  also  arise  in  industrial  accidents.  The  ground  and  dynamic 
impulse  generated  by  such  explosions  is  a main  cause  of  damage; 
hence,  it  is  important  to  determine  the  influence  of  the  fuel  properties, 
the  cloud  geometry,  and  the  size  of  the  initiating  charge  upon  the 
impulse  which  is  generated.  The  analytical  determination  of  the 
impulse  generated  by  a fuel-air  explosion  forms  the  subject  of  the 
present  paper. 

The  shape  of  the  fuel  cloud  and  the  fuel  distribution,  which  will 
have  important  effects  upon  the  ground  impulse,  will  generally  be 
quite  complex,  depending  on  the  method  of  fuel  dispersal.  Analytical 
determination  of  ground  impulse  therefore  is  possible  only  if  simpli- 
fied models  of  fuel-air  explosions  are  used.  Use  of  appropriately 
formulated  models  makes  it  possible  to  rapidly  estimate  the  effects 
of  the  fuel  properties  and  cloud  geometry  upon  the  impulse  and  also 
provides  a basis  of  assessing  experimental  data  from  fuel-cloud 
explosions. 
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Hence,  it  is  assumed  that  fuel  is  distributed  uniformly  through 
the  cloud  with  Q,  the  combustion  energy  released  per  unit  mass  of 
fuel-oxidizer  mixture.  Clouds  with  planar,  cylindrical  and  spherical 
symmetry  are  considered.  It  is  known'1'2,3,4'  that  the  initiating 
blast  wave  Is  dominant  at  first  but  that  later  the  explosion  consists 
of  a Chapman-Jouguet  detonation  propagating  through  the  fuel-cloud. 
The  transition  between  blast  wave  and  detonation  behavior  occurs 
near  a critical  radius  r,  within  which  the  total  combustion  energy 
release  equals  the  energy  E„  of  the  initiating  blast  wave.  This 
condition  leads  to  the  relation 

(1> 

where  v = 1,  2,  3ando„=2,  2v,  4 v for  planar,  cylindrical,  and 
spherical  symmetry.  The  flow  in  the  transition  region  when  the 
explosion  radius  r,  is  of  the  order  of  r.  Is  very  complicated. 
However,  it  is  shown  by  Nicholls  et  al.'  ' that  a zeroth  order 
model  in  which  self-similar  blast  wave  and  detonation  wave  theory 
are  used  for  r9  < r,andrB  > respectively,  yields  analytical 
results  in  reasonable  agreement  with  experiments  carried  out  in  a 
special  pie-shaped  shock  tube,  designed  to  simulate  an  explosion 

with  cylindrical  symmetry. 

This  zeroth  order  model  also  is  the  basis  of  the  Impulse  cal- 
„„iatinna  described  below,  and  leads  to  results  which  permit  rapid 


computation  of  the  impulse  using  universal  impulse  functions.  In 
the  planar  and  cylindrical  cases  it  is  at  first  assumed  that  the  fuel- 
cloud  extends  from  the  ground  plane  to  an  infinite  height  or,  what  is 
equivalent,  that  the  explosion  is  confined  by  a non-yielding  plane  which 
is  parallel  to  the  ground.  With  this  model  the  influence  of  the  fuel 
used,  the  fuel-oxidizer  mixture  ratio,  and  the  blast  energy  upon  the 
impulse  generated  by  the  wave  is  readily  determined.  The  total 
ground  impulse  and  the  static  dynamic  impulse  at  a fixed  distance 
from  blast  center  generated  by  the  explosions  of  such  confined  clouds 
are  considered  in  Parts  B through  G. 

Actually,  of  course,  the  upper  surface  of  the  fuel-clouds  is 
bounded  by  inert  air  or  oxidizer  rather  than  by  a non-yielding  surface. 
Thus,  as  the  detonation  propagates  through  the  cloud  an  expansion 
wave  will  eat  into  the  high  pressure  combustion  products  behind  the 
wave.  This  side  relief  will  cause  considerable  reduction  in  the 
ground  impulse  generated  by  the  wave  and  will  depend,  mainly,  on 
the  height  of  the  fuel*cloud.  Simplified  models  for  estimating  this 
effect  are  considered  in  Part  H. 

B.  THE  TOTAL  GROUND  IMPULSE 

The  pressure  p,  within  the  leading  shock  of  a FAE  (fuel-air 

explosion),  i.  e.,  in  the  region  r < r , will  be  a function  of  the  radial 

s 

distance  r from  the  blast  center.  If  p^  is  the  ambient  pressure  ahead 


of  the  wave,  then  the  ground  impulse  generated  by  the  FAE  per 
unit  area  at  a given  point  during  the  time  interval  dt  is  (p  - Pj)dt. 
The  total  impulse  dl  generated  in  the  region  r < r during  dt  is 


then  given  by 


( rs(t) 

* j [pfr.t)  - P^J 


<jv  rV  * dri  dt 


The  index  v = 1 for  planar  waves  but  v = 2 for  both  cylindrical  and 
spherical  waves  since  each  exerts  pressure  on  a circular  area. 
Fuel-air  explosions  with  planar,  cylindrical  and  spherical  symmetry 
are  shown  in  Figure  1.  The  tcAal  ground  impulse  I (t)  at  time  t after 
initiation  of  the  secondary  blast  will  then  be 


I(t)  = o 


i 

I [p(r,t)  - pj 

i 0 


v-1  , .. 

r dr  dt 


It  is  important  to  recognize  that  in  the  planar  cas&I(t)  is  the  impulse 
applied  to  a strip  of  unit  width. 

In  general,  evaluation  of  the  impulse  from  Equation  (3)  will 
require  extensive  numerical  integration.  The  integral  (3)  can, 
however,  be  greatly  simplified  for  the  self -similar  blast  waves 
and  detonation  waves  which  are  used  in  the  zeroth  order  FAE  model 
described  in  Part  A. 
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C.  BLAST  WAVE  IMPULSE 

When  r < r*  the  FAE  is  thus  replaced  by  a strong  blast  wave. 
The  self-similar  solution  for  a strong  blast  wave  is  discussed  in 
detail  by  Sedov (5)  and  Taylor (6).  The  key  assumptions  are  that  the 
shock  pressure  ratio  (P2/  Pj)  » 1,  that  the  release  of  explosive 
energy,  EQ,  is  instantaneous,  and  that  the  fluid  is  a perfect  gas, 
with  constant  specific  heats.  Then  it  can  be  shown  that  the  shock 
radius  r , and  the  fluid  velocity  v2,  pressure  p2,  and  temperature 

T-  immediately  behind  the  shock  front  are  given  by 

21 


V2=F 


P*  = 


T2  = 


4 

+ + 1)  VPij 

8 pi  Kf 

* 2)2  (r  ♦ l)  W 

8(y-t)  /M 

■ + 2)2  (r  + l)z  \°pl/ 


-v/ 1/+2 


(4) 


-2v/  v+2 


-2v/  v+2 


E is  the  blast  energy  released  instantaneously  per  unit  area,  length, 
o 

or  at  a point  for  v = 1,  2,  3.  The  function  a (y,  v),  which  is  derived 

from  the  requirement  that  the  total  energy  for  r < rg  must  equal 

E is  shown  in  Figure  2.  The  dimensionless  pressure 
o 


(5) 


l 


Pbfc,  v)  - p/ px  /^=  1 1 

b 1 \t2/  \ 2/  (v  + 2)2  (y  + 1)  XZ 


is  a function  of  only  the  normalized  radius  X = r/ r , the  geometric 

s 


in^ex  v,  and  the  ratio  of  specific  heats  y. 

Neglecting  p^,  and  introducting  Equations  (4)  and  (5)  in 
Equations  (3)  reduces  the  impulse  integral  to  the  form 


with 


(6) 


0(r 


1 

- ^ = f pb 


v ) XV+*  dX 


The  integral  /3(y,  v)  depends  only  on  the  rati.)  of  specific  heats  y and 
the  geometric  index  v,  and  can  be  evaluated  from  the  self-similar 
blast  wave  solution  which  is  presented  in  analytical  form  by  Sedov 
For  spherical  blast  waves,  which  act  on  a circular  area,  the  total 
impulse  is  given  by 


(5) 


yt)  = -j7T  Pl{Eo/aPl)^5  0(y,3)  t3^5 


with 


CO 


0(y 


1 

,3)=  j 


Pb(X,3)A  dA 


I 1 


The  functions  /3(y,  u)  for  v = 1,2,  3,  are  plotted  in  Figure  3.  The 
impulse  from  planar  and  cylindrical  waves  increases  more  rapidly 
with  time  than  for  a spherical  wave.  This  result  is  to  be  expected 
since  the  spherical  geometry  provides  more  pressure  relief  than 
the  confined  planar  and  cylindrical  geometry.  ' 

Equations  (6)  and  (7)  also  suggest  the  introduction  of  a nor- 
malized impulse  I defined  by 


= i ''■M 


•b  ■ yti/yvp/75 13/6  - 


(8) 


v = 3 


The  variation  of  the  normalized  impulse  I with  the  ratio  of  specific 
heats  y is  shown  in  Figure  4.  Using  Equation  (4)  the  time  t in  the 
expressions  for  I(t)  can  be  replaced  by  the  shock  radius  so  that 


v+2 

2 ; ‘'  = 1-2 

(9) 

yrs)  = T"(E</Va)1/2^s3/2  ! 9*3 

D.  DETONATION  WAVE  IMPULSE 


When  rg  > r*,  the  self-similar  solution  for  the  tlow  behind 
a Chapman -Jouguet  (CJ)  detonation^  is  used  to  represent  the  FAE. 
The  total  ground  impulse  generated  by  a CJ  detonation  is  also  given 


H mm 


by  Equation  (3);  however,  the  ambient  pressure  is  no  longer 

negligible  compared  to  the  pressure  Pg  behind  the  detonation  front. 

As  before,  a dimensionless  pressure  P , defined  by 

d 


Pd(*»  v)  = p/px 


(10) 


is  introduced.  In  contrast  to  the  strong  blast  wave,  a detonation 
front  propagates  with  a constant  velocity  C so  that 


rB  " Ct  (11) 

Precise  determination  of  the  pressure  ratio,  temperature 
ratio,  density  ratio  and  propagation  speed  of  CJ  detonations  re- 
quires the  determination  of  the  equilibrium  composition  of  the 
products  of  combustion,  and  has  been  carried  out  numerically  using 
the  computer  code  of  Gordon  and  McBride^.  Detonation  param- 
eters for  methane-air  and  MAPP  gas-air  mixtures  computed  with 
this  code  are  shown  in  Figures  5 and  6 and  are  typical  of  hydro- 
carbon-air detonations.  Once  the  propagation  velocity  C,  detonation 
Mach  number  the  ratios  of  specific  heats  and  y^,  and  the 
molecular  weight  Wjf  of  the  combustion  products  have  been  deter- 
mined numerically,  other  detonation  properties  can  be  determined 
from  the  perfect  gas  Hugoniot  Equations  with  remarkable  accuracy. 

The  equations  for  mass,  momentum,  and  energy  conservation  across 

13 
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/c  o\ 

the  detonation'  ’ ',  the  perfect  gas  equation  of  state  for  the  fuel- 
air  mixture  ahead  of  the  wave  and  the  combustion  products  behind 
the  wave,  and  the  CJ  condition  that  the  downstream  velocity  is  sonic 
relative  to  the  wave  then  lead  to  the  relations: 


fa. 

Pi  Vo 


f * -^y1 

V nMd/ 


P2  l+V1Md 


pl  *2  * 1 


> ■ 4,!  ■ ■) 


T 


Tj  2rlr2  ® 


<Vl)  A ^2 


(r i"1)  (y2+l) 


with 


(12a) 

(12b) 

(12c) 

(12d) 

(12e) 


(yri)Q 


(yri)Q 


The  sign  = denotes  equations  which  are  valid  only  when  ^ » 1,  a 
condition  which  is  valid  for  all  detonations  under  consideration  here. 
For  example,  given  C and  y ^ from  the  numerical  computations, 
Equation  (12d)  can  be  used  to  compute  the  heat  release  Q, 
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which  also  occurs  in  the  expression  (1)  for  critical  radius  r*. 


Introduction  of  Equations  (10),  (11),  and  the  similarity  vari- 
able X = r/rg  into  the  equation  for  the  total  impulse  yields  the  result 


Vi c 


4 3 

2tt  p1  C t 

V‘> 3 — 
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: 
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where 


Vv)  = ( Pd 


0(r2,3)=  J Pd(X,3)\3dX 


v = 1,2 


v = 3 


v)\V+1dK  ; v - \f2 


; v = 3 


Using  Equation  (11)  the  impulse  Id  can  also  be  expressed  in  terms 


of  the  detonation  radius  r with  the  result 

s 


r,  . %pic 

W" rl 


, v+i  r 1 

TT — ' v 

L v*i  Md  J 


= 1,2 


2ir  p1  C: 

Id(r»)  * — r- 


v = 3 
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The  variation  of  Pd(X,  i/),  and  hence  the  self -similar  solution 
for  the  flow  behind  a CJ  detonation  must  be  known  to  compute  the 
integrals  6(y2>  v).  This  self -similar  solution  cannot  be  found  in 
analytical  form,  and  as  indicated  by  Sedov ^ the  key  step  is  inte- 
gration of  the  non-linear  differential  equation 

dz  _ z[2(V  - l)2  + ( v - 1 ){y2  - 1)  V(V  - 1)  - 2z] 

57  (16) 

V[  (V  - l)2  - vz] 

The  dimensionless  velocity  V(X)  and  the  variables  z(X)  are  defined 
by 

V Pdfc> 

v(x)  = (FTtT  ; z=y2W^r)  > = ^ (17) 

Integration  of  Equation  (16)  must  be  carried  out  numerically  and 
starts  at  the  point 


immediately  downstream  of  the  detonation  front.  The  solution  z(V) 

thus  depends  on  both  and  the  detonation  Mach  number  Md;  however, 
2 

since  Md  » 1,  this  latter  dependence  is  very  weak.  The  detona- 
tion front  is  followed  by  an  isentropic  expansion  through  which  the 
velocity  of  the  combustion  products  decreases  to  zero  at  the  boundary 
of  a stationary  core  region  corresponding  to  the  point  V = 0,  z = 1,  a 
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singular  point  of  Equation  (16).  Details  of  the  numerical  solution  of 
Equation  (16)  are  described  in  Reference  9.  All  other  parameters 
of  the  flow  behind  the  detonation  are  readily  determined  once  the 
variation  of  z with  V is  known(5,9).  Typical  velocity,  pressure,  and 
temperature  profiles  behind  a CJ  detonation  are  shown  in  Figure  7. 

The  function  6(y2,  v)  has  been  calculated  for  both  air -methane 
and  air-MAPP  mixtures.  For  the  two  fuels  considered,  y2  ranges 
between  1. 16  and  1. 31  while  Md  lies  between  4. 0 and  5. 0,  and  over 
this  range  5(y0,  v)  is  almost  independent  of  y2  and  the  fuel-oxidizer 
mixture  ratio.  In  fact,  6(y2,  v)  depends  only  on  v for  methane  and 
MAPP,  and  the  average  value  of  6(y2,  v)  is  given  in  Table  I below. 

The  values  of  y2  and  Md  for  MAPP  and  methane  are  typical  of  most 
hydrocarbon-air  detonations;  hence,  the  values  of  6(y2,  v)  in  Table  I 
can  be  used  to  compute  the  impulse  for  hydrocarbon -air  CJ  detonations 

in  general  with  reasonable  accuracy. 

TABLE  I.  AVERAGE  VALUE  OF  «(r2>  v) 

I 

v 6 ^2*  ^avg 

1 0. 190 

2 0. 138 

3 0. 135 


The  quantities  Fj(t)/ 1^+1  for  v = 1, 2,  and  Tj(t)/ 1^  for  v = 3 
depend  on  p^,  C,  and  (Equation  13)  and  so  vary  appreciably  with 
the  type  of  fuel  and  with  the  equivalence  ratio  0.  With  these  impulse 
functions,  which  are  plotted  for  methane-air  and  MAPP-air  mixtures 
in  Figures  8 and  9,  Id(t)  is  readily  determined  as  a function  of  0.  The 
functions  Id^)/tt’+1,  Id(t)/t^,  also  provide  a basis  for  determining 
the  influence  of  the  fuel  composition  and  mixture  ratio  upon  the  total 
ground  impulse.  From  Figures  8 and  9,  for  instance,  it  can  be  seen 
that  MAPP-air  mixtures  provide  a somewhat  higher  impulse  than 
methane-air.  It  is  also  of  interest  that  the  peak  impulse  occurs 
for  rich  mixtures  with  0 = 1. 2 and  0 = 1. 5 for  methane-air  and 
MAPP-air,  respectively. 

E.  GROUND  IMPULSE  OF  FAE 

Both  the  blast  wave  and  the  CJ  detonation  results  must  be 
used  to  determine  the  total  ground  impulse  f of  a FAE.  When 
rs  < r*,  F is  given  by 

•<\>  -~W  <19> 

for  the  zeroth  order  FAE  model  described  in  the  Introduction.  The 
blast  wave  and  CJ  solutions  are  patched  together  when  r = r4. 

This  means  that  the  entire  flow  field  corresponds  to  that  of  a blast 
wave  when  r < r#,  and  to  that  of  a CJ  detonation  when  r > r . . 

a S “ 

Accordingly  the  ground  impulse  is  approximated  by 
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I(rs)=Id(rs)  “ ld{r*]  + Vr*>  (20) 

Thus  when  rg  = r+,  T(rg)  equals  the  blast  wave  ground  impulse,  while 
with  rg  > r+,  T(rg)  equals  the  ground  impulse  generated  by  the  blast 
wave  up  to  rg  = r^,  plus  the  ground  impulse  generated  by  the  CJ 
detonation  for  rg  > r*.  T^r^  and7d(rg)  are  given  by  Equations  (9) 
and  (15),respectively. 

As  an  illustration  the  total  ground  impulse  will  be  determined 
for  a cylindrical  FAE  through  an  air-MAPP  mixture  at  atmospheric 
pressure  with  equivalence  ratio  0 = 0.  563.  For  this  mixture: 


Md  = 4’69  r2=1*23 


C = 5320  ft/  sec 


y j * 1.  30  Pj  = 2. 29  x 10  3 slug/  ft3  P1  = 14.  7 psia 

Experiments  ^ indicated  that  the  critical  blast  energy,  E , needed 
to  detonate  this  mixture  was  4.  38  x 105  (ft-lb)/ft,  equivalent  to  150 
grams  of  Detasheet  per  foot.  For  these  conditions 


2Eo  - »> 

’yl»lMd2 


= 1.54  ft 


where  Equation  (12c)  has  been  used  to  express  Q in  terms  of  M 

d' 

With  r < r. 


P(yv  2)  rgZ  = 9.  36  rsZ(]b-sec); 


1<rs)  = W = 2” 


[^ij 


r in  ft 
s 


or  in  terms  of  time 


I(t)  = yt)  = 1. 17  x 105  t (lb-sec) 


When  r > r* 
s * 


_ _ _ _ 27 r p«  C 

■ w -1d(r*)  +\(r«)  — r- 


L *lMd  J 


, 3 
(rs  -3 


1/2 


/VlY  2 3 

2*  2)  r/  = (3.  06  r^  + 11.  l)(lb-sec) 


as  a function  of  time  the  impulse  becomes 


I(t)  = 


2tt  pj  C 


M 


(t3  -t.3)  + 2jt  E 0/ a)t* 


= (4. 60  x 10U  t3  + 19.-3)  (lb-sec) 


Here  t+  corresponds  to  the  time  at  which  the  blast  wave  radius 
reaches  the  critical  value  r^  so  that 

1/2 


-4 


For  the  example  above  t*  = 1. 88  x 10  sec. 


J 


» 
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The  total  FAE  ground  impulse  is  plotted  as  a function  of  both 
rg  and  t in  Figures  10  and  11  for  the  illustrative  example  computed 
above.  It  is  of  interest  to  note  the  small  value  of  r*  in  the  case 
treated  above,  indicating  that  a fuel-cloud  with  a radius  of  the  order  of 
10  feet  would,  in  this  case,  be  dominated  by  the  CJ  detonation. 

This  result  also  suggests  that  the  blast  wave  phase  of  the  FAE  can 
in  some  cases  be  neglected  entirely. 

F.  GROUND  IMPULSE  AT  A FIXED  DISTANCE  FROM  BLAST 
CENTER 

In  addition  to  the  total  impulse  it  is  of  considerable  interest 
to  determine  the  ground  impulse  per  unit  area  which  is  generated 
by  a FAE  at  some  fixed  distance  r from  the  blast  center.  If  r > r* 
the  impulse  is  generated  entirely  by  the  passage  of  a CJ  detonation 
if  the  zeroth  order  model  of  the  FAE  is  used,  and  the  present  dis- 
cussion will  be  restricted  to  this  case. 

The  impulse  ^ per  unit  area  at  distance  r from  blast  center 
is  given  by  the  integral 

t 

Jl(r,t)=  j £p(r,t)  - p^j  dt  (22) 

s 

The  lower  limit,  t , of  this  integral  is  the  time  at  which  the  wave 

s 

front  first  arrives  at  r and 
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Ground  Impulse 


*s  - r/C 

Obviously  <1.  0 when  t < t . In  Equation  (22)  and  in  the  analysis 

s 

below,  r plays  the  role  of  a parameter  and  the  time  t is  the  only 

variable.  From  Equation  (10)  it  follows  that  (p/pg)  is  a function 

of  the  similarity  variable  X = r/ r = r/Ct.  Then  introducing  X as 

s 

the  variable  of  integration,  Equation  (22)  can  be  expressed  in  the 
form 

o,> 

The  integral  in  (23)  depends  on  and  weakly  upon  the  detonation 
Mach  number  M^.  Since  y ^ in  any  case  lies  in  a very  narrow 
range  this  integral  can  essentially  be  treated  as  a universal  im- 
pulse function.  Thus  it  is  useful  to  define  a function  U(X)  by 

W|-i m 

The  function  U can  also  be  considered  a function  of  a dimensionless 

time  r = (t/t  ) = x’1  so  that 
s 

T pi 

u(t)  = + zr  (T  - (24a) 

P2s  P2 

As  indicated  in  Part  D above  the  detonation  is  followed  by  an 
expansion  wave  in  which  the  fluid  velocity  drops  to  zero  at  some 
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\ = \ The  expansion  is  followed  by  a stagnant  core  from  X = X2 
to  blast  center  X = 0 where  the  pressure  is  constant  and  equal  to 
some  value  Pg.  Surprisingly  (pg/  p2)  0.  35  and  Xg  ~ 0.  5 for  v = 1, 2, 3 

and  for  the  range  of  y 2 typical  of  hydrocarbon  fuels.  Taking  this 
behavior  into  account  it  is  possible  to  express  U(X)  in  a somewhat 
simpler  form,  especially  for  X < Xg.  Thus  when  X2  K ^ < 1.0 


U(X)  = 


J?© 


dX 


(25) 


When  0 < X < X, 


(26) 


or  expressed  in  terms  of  t with  r 2 = Xg 


-1 


U(t)  = U(t2)  + — (t  - r2) 


(27) 


Since  U(t2)  - 0(1/  2)  it  is  clear  that  for  large  times  r » t^,  U(t) 
will  be  dominated  by  the  second  term  on  the  right-hand  side  of 
Equation  (27). 

It  is  useful  and  revealing  to  write  down  the  expressions  for  the 
impulse  in  terms  of  the  universal  impulse  function  U(t).  Thus 


= 0 


T < 1 


(28a) 
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T > 1 


(28b) 
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and  when  r > r 2,  Equation  (28b)  can  also  be  written  in  the  form 


P3.  . pl 


v2=p2§u<T2)+e(T-T2>-i:<T-M 


(29) 


-2  “ -2 

Generally  U(t)  » (Pi/Po)(T  - 1)  and  then  it  follows  from  Equation 
(28b)  and  Equations  (12)  that 


- Md_r  Pl 
'2 


J * tZtv  uw 


(30) 


and  Equation  (30)  shows  in  a very  simple  way,  how  *J)  depends  on  r, 
the  distance  from  blast  center,  and  the  detonation  parameters  Md,  pj, 
and  y^. 

There  remains  the  problem  of  evaluating  the  universal  function 
U(t).  For  plane  waves  with  v = 1,  (p/p2)  can  be  expressed  in  analy- 


tical form^’  ^ and  for  X2  < X < 1 is  given  by 


JL  = 

P2 


I c ^2  ” ^ 

= (x  - 1) 

a2y2  + 1 ] 


^2y2/(y2-l) 


y2  - 1 

= 1+-V-&-D 

L y2 


Jy2/(y2-V 


(31) 


The  second  form  in  Equation  (31)  follows  from  using  (p2/  p^)  = (y2+l)/y2 

33 


, . . ...  _ 


j 


. ... — — i. 


from  Equation  (12a).  With  v = 1 it  follows  that  X = (1/2) 


and 


-2y2/"/2'1 


p3  y2  + 

P^=  ~^~l  (32) 

With  Equations  (31)  and  (32)  u(r)  is  then  readily  determined.  For 
cylindrical  and  spherical  CJ  detonations  with  v - 2,  3 it  is  necessary 
to  use  numerically  computed  values  of  (p/p2)  as  shown  Figure  , 


The  universal  impulse  function  U(t)  computed  for  v = 1 
in  Figure  12. 


> 2,  3 is  shown 


G'  blastMcenterLSE  AT  A FKED  distance  ™om 

Any  body  or  obstacle  in  the  path  of  the  FAE  will  be  subjected 

to  a dynamic  force  due  to  the  leading  shock  front  and  the  region  of 

high  velocity  immediately  behind  the  shock.  As  indicated  by  Erode 
(11)  * 
ct  al  the  blast  load  will  consist  of  the  forces  generated  when  the 

shock  strikes  the  obstacle,  reflects  and  eventually  flows  past  the 

structure.  After  sufficient  time  has  elapsed,  it  seems  reasonable 

to  assume  that  the  flow  will  be  quasi-steady,  i.  e,  that  the  flow  behaves 

like  a steady  flow  with  a free  stream  velocity  equal  to  v,  the  velocity 
behind  the  shock  or  blast  wave. 

During  the  quasi-steady  portion  of  the  flow  the  force  on  the 
body  due  to  the  induced  air  flow  will  be 
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Universal  Impulse  Function 


F =£La  c 
d 2 Af  CD 


where  Af  is  the  frontal  area  of  the  body,  and  CD  is  the  drag  coefficient. 
The  dynamic  impulse  F applied  to  a body  at  distance  r from  blast 


center  will  then  be 


Fddt  = 


u «— * 

)V  A 

T Af  °D  dt  = CD  Af 


where 


4 - [ 

v I 2 


and  CD  is  an  average  value  of  the  drag  coefficient.  The  evaluation  of 
JIv,  the  dynamic  impulse  per  unit  area  and  drag  coefficient  in 
terms  of  universal  functions  will  be  considered  below.  Once  Af  and 
CD  are  known  is  then  readily  determined  from  In  the  absence 

of  more  specific  information  the  assumption  that  CD  = 1 can  often  be 
used  successfully  for  preliminary  results. 

It  must  be  emphasized  thatlv,  given  by  Equation  (34\does  not 
take  into  account  the  unsteady  phenomena  which  occur  when  the 
leading  shock  first  reaches  the  obstacle.  At  the  first  instance  the 
shock  is  reflected  and  the  pressure  at  the  obstacle  surface  reaches 
the  reflected  shock  pressure  which  may  be  much  greater  than  the  C J 
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pressure,  p2-  If  L is  the  characteristic  size  of  the  obstacle  then 

it  appears  reasonable  that  a characteristic  time  r for  the  initial 

transient  period  will  be  of  the  order  of  (L/  C).  Since  the  decay  of  the 

induced  velocity  v behind  a CJ  detonation  occurs  in  the  range 

(1/ 2)  < X < 1 the  time  r during  which  the  dynamic  force  F . 

i d 

acts  on  the  obstacle  will  be  of  the  order  r/C.  Thus  the  ratio  of 
(Tt/  Td)  = (L/ r)  will  be  small  provided  (L/r)  « 1,  i.  e.,  provided 
the  obstacle  size  is  small  compared  to  the  distance  from  blast 


center. 


Equation  (35)  for  can  also  be  expressed  in  terms  of  the 


similarity  variable  X with  the  result 


X2  < X < 1 or  tg  < t < t2 

and  V^v  remains  constant  at  the  value  corresponding  to  X_  when 

x < X2.  Using  the  Hugoniot  Conditions  [Equation  (123  with 
2 

ylMd  ^ ^ makes  it  possible  to  express  Equation  (36)  in  the 
following  more  useful  form: 


37 


and  Uv(X)  or  Uy(T)  is  again  a universal  impulse  function  which  is 
essentially  independent  of  detonation  properties.  The  universal 
impulse  function  Uv(A)  is  shown  for  v = 1,  2,  3 in  Figure  13. 

H.  THE  INFLUENCE  OF  SIDE  RELIEF 

The  analysis  of  impulse  described  above  considers  cylindrical 

or  plane  waves  that  are  of  infinite  height  or,  what  is  equivalent, 

that  are  confined  by  a non-yielding  surface  parallel  to  the  ground. 

Actually  such  clouds  are  bounded  by  inert  air  or  oxidizer.  As  the 

detonation  propagates  through  the  cloud  an  expansion  wave  will  thus 

pr  opagate  through  the  combustion  products  causing  a considerable 

reduction  in  the  impulse  which  is  generated. 

The  influence  of  side  relief  on  a planar  detonation  has  been 

considered  by  Sichel(12),  Dabora.et  al(13',  and  Sommers  and 
(14) 

Morrison  , and  this  geometry  will  also  be  considered  hero. 
Analysis  of  a planar  wave  with  side  relief  makes  it  possible  to 
establish  the  influence  of  cloud  thickness  upon  impulse.  Cylindrical 
waves  witc  side  relief  are  more  complex  to  analyze;  however,  for 
large  values  of  the  wave  radius  the  flow  near  the  wave  front  will 
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approach  that  behind  a plane  wave.  In  the  analysis  below  it  will 
be  assumed  that  the  initiation  of  the  detonation  wave  is  almost  in- 


stantaneous or  that  r*  « r . The  sample  computations  in  Part  E 

¥ s 

indicate  that  this  is  not  an  unreasonable  assumption. 

Propagation  through  a cloud  of  height  h will  be  considered  as 
shown  in  Figure  14.  Initially,  side  relief  will  not  influence  the 
variation  of  the  pressure  on  the  ground  plane.  During  an  initial 
period  of  duration  t^  the  flow  behind  the  detonation  changes  from  that 
of  the  self-similar  solution  described  in  Part  D to  the  quasi-steady 
flow  shown  in  Figure  15.  The  duration  of  this  initial  period  will  be 
of  the  same  order  as  the  time  required  for  a sound  wave  to  travel 
the  width  of  the  fuel-cloud,  i.e.,th  ~ 0(h/a2)  where  ag  is  the  speed 
of  sound  immediately  behind  the  detonation  front.  To  simplify  the 
analysis,  the  self-similar  CJ  wave  pressure  distribution  will  be 
used  to  compute  the  wave  impulse  during  the  initial  period  t < t^. 
During  time  t^,  the  detonation  will  travel  a distance  L from  blast 
center  given  by 

L = CL=cf  (38) 

II  ^ n 


From  the  CJ  conditions  it  follows  that  (C/ a2)  = (p2/  p^)  so  that 


L P2  J2  + 1 

h"^r  ^2 
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(39) 
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Since  = 1. 2 for  the  products  of  combustion,  it  follows  that  the 
influence  of  side  relief  must  be  considered  after  the  detonation 
has  travelled  a distance  equal  to  about  twice  the  height  of  the  fuel- 
cloud. 

When  t > the  flow  behind  the  wave  will  be  as  pictured  in 


Figure  15,  which  shows  the  CJ  detonation  in  a coordinate  system 
fixed  to  the  wave.  The  CJ  detonation  is  followed  by  an  expansion 


wave  across  which  the  pressure  drops  from  the  CJ  pressure  P to 

e2 


PCg  while  the  combustion  products  turn  through  an  angle  6.  The 
molecular  weight  of  the  inert  ^ ^ will  in  the  case  of  air  be  of  the 
same  order  as  ^ j,  the  molecular  weight  of  the  unburned  ex- 
plosive. The  CJ  detonation  generally  induces  an  oblique  shock 

,(12) 


wave  in  the  inert  gasv  ' across  which  the  flow  is  deflected  through 
the  angle  6.  The  interface  conditions  are  determined  by  the  re- 
quirements that  the  oblique  shock  wave  and  the  expansion  wave 
both  turn  the  flow  through  the  same  angle  6,  and  that  the  pressures 
Pj2  and  PGg  behind  the  shock  and  the  expansion  be  the  same.  The 
expansion  wave  is  reflected  from  the  ground  plane  where  the  pres- 


sure drops  from  the  CJ  value  Pp„  to  Pp  at  a distance  f from  the 

e2  e4 


detonation  front.  The  region  r < (r  - £ ) is  neglected  for  the  pur- 

s 


pose  of  computing  the  ground  impulse.  Here,  the  pressure  will 
generally  be  only  slightly  above  atmospheric  and  will  ultimately  be 


1 


4 

1 1 
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relieved  by  further  reflections  from  the  interface  between  the  com- 


bustion products  and  the  inert  gas.  The  interaction  at  the  inert- 
explosive  interface  causes  a reduction  in  the  propagation  velocity 
C which  is  proportional  to  the  ratio  of  the  reaction  zone  thickness 
A to  the  wave  height  h.  This  effect  will  be  neglected  here. 

Since  = 1.0,  the  deflection  angle  6 and  the  Mach  Number 


downstream  of  the  expansion  wave  are  related  by 


(15) 


6 = 


ye2  + 1 


*|1  / 2 


ye2  “ 1 


tan 


i 1/2 


1/2 


(40) 


From  the  oblique  shock  relations 


6 = tan 


-1 


2 cot  9 


M.-2  sin2  0-1 
ll 


(41) 


Mil  ^il  + COS  2^  + 2J 

and  by  equating  Equations  (40)  and  (41)  one  relation  between  and 
9,  the  shock  angle  (Figure  15)  is  obtained.  M.  ^ and  can  be 
determined  from  the  detonation  velocity  and  the  properties  of  the 
inert  gas.  The  pressure  ratio  across  the  oblique  shock  will  be 


P. 


i2 


= 1 + 


(Mil2  si”2  e - x) 


(42) 


The  flow  across  the  expansion  wave  is  isentropic.  Then  using 
Equation  (12b)  for  the  pressure  ratio  across  the  detonation  it  is 


readily  shown  that 
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Since  OP^/^l  = Equations  (42)  and  (43)  provide  the  addi- 


tional equation  required  to  determine  M „ and  0.  Once  M « is 

eo  eo 

determined,  6,  Pgg  and  the  details  of  the  expansion  wave  are  readily 

established.  The  pressure  variation  along  the  ground  plane  due  to 

the  reflected  expansion  wave  can  be  determined  using  the  method  of 

characteristics  as  described,  for  instance,  in  Reference  15. 

Following  the  discussion  above  the  impulse  I(r  ) will  be  given 

s 

by 

f(rs)  = W ; rs  < L = ah  <44> 

withTd(rg)  from  Equation  (15).  There  will  be  cases  in  which  the 
reflection  length  SL  > L.  In  this  case  confmed  wave  impulse  1.(r  ) 

u S 

will  be  used  until  the  wave  has  moved  distance  l from  blast  center, 
i.  e. , until  the  expansion  wave  reflection  pattern  is  completely 
established.  With  rg  > L or  f the  total  ground  impulse,  Ig,  of  a 

planar  detonation  with  side  relief  will  be 
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(P  -P^dr  (t  -th)  + Id(th) 


(45) 


The  integral  in  Equation  (45)  is  more  conveniently  expressed  in 

* 

terms  of  the  variable 


where  x is  the  distance  of  a given  point  from  the  wave  front  as  shown 


in  Figure  15.  Then  Equation  (45)  can  be  written  in  the  form 


Since  h is  the  only  characteristic  length  to  enter  the  interactions 

between  the  combustion  products  and  the  inert  boundary,  h can  be 

used  to  scale  the  dimensions  of  the  interaction  region.  It  follows 

that  the  wall  pressure  variations  P/  P « will  depend  only  on  the 

e* 

dimensionless  variable  x/h  for  a given  fuel-oxidizer  mixture.  A 
consideration  of  the  expansion  wave  behind  the  detonation  shows  that 


Herice,  the  integral  in  Equation  (46)  which  will  be  denoted  by  U ($) 

s 

so  that 
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mem 


h r 

V*>  -J  |j 


P r 1 

/-(U-p1  d* 

e2  e2 


will  depend  only  on  the  equivalence  ratio  $ for  a given  fuel. 

It  is,  perhaps,  most  convenient  to  write  the  expression  for 


I in  the  form 
s 


ye2  + 1 

r > — h 

6 ye2 


(46a) 


The  fuel  properties  will  determine  P g/ Pj,  C,  and  Ug($).  The 

height  h enters  the  expression  for  I in  two  ways.  First  h appears 

s 

explicitly  in  the  coefficient  of  U ($),  then  h determines  the  distance 

s 

L the  detonation  travels  before  side  relief  becomes  important.  If  R 
is  the  outer  radius  of  the  fuel-cloud  and  if  h > y^R/  (ye2  + *)>  the 
foregoing  considerations  indicate  that  side  relief  will  not  affect  the 
total  ground  impulse. 

The  function  U ($)  has  been  computed  for  MAPP-air,  and 
s 

methane-air  mixtures.  Over  the  range  0. 6 < $ < 1.6  the 
Prandtl  Meyer  angle  6 and  the  shock  angle  6 were  found  to  remain 
essentially  constant  at  6 = 20°,  6 = 29°  for  both  mixtures.  As  a 
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consequence  the  ratio  H/h  = 2.  90  and  is  also  constant  so  that  the 
integral 


(4  ) d£  = 1.  3885 


and  is  also  constant  for  the  above  two  fuels.  The  impulse  function 

U ($)  then  has  the  remarkably  simple  form 
s 


U ($)  = 1.3885  - f — 

S \ J pe2 

\ 0 


* P!  ^ — pt"" 

= 1.3885  - £ — = 1.3885  - 2.90—  (48) 

hpe2  Pe2 

It  must  be  emphasized  that  this  result  is  valid  only  for  CH^-air  and 

MAPP-air  mixtures  in  the  range  0. 6 < $<  1.6.  However,  it  is 

likely  that  Equation  (48)  will  remain  valid  for  other  hydrocarbon -air 

mixtures  over  this  range  of  equivalence  ratio. 

For  large  values  of  r /h  the  total  impulse  I generated  by  a 

s s 

planar  wave  with  side  relief  is  approximately 


P,  hr  /P  _\ 

- - V(l $ Us<*> 


while  the  total  impulse  of  a planar  confined  CJ  detonation  will, 
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from  Equation  (15)  be  approximately 

V^e2  + i)(^)[a<ye2,  i)--^  <*» 

— _ 2 — 

Thus,  I oc hr  and  I.  ccr  so  that  I.  increases  more  rapidly  with  r 
ssds  d s 

than  I . This  result  is  to  be  expected  since  a core  of  high  pressure 

stagnant  gas  is  maintained  behind  the  confined  detonation  front;  whereas 

this  high  pressure  region  is  destroyed  when  there  is  side  relief.  Thus 

when  r / h » 1 
s' 


I 

■ 


L 


I /L  oc 
s d r 


(51) 


s 


When  (£/ r ) « 1,  the  flow  behind  a cylindrical  detonation 

D 

with  side  relief  will  approach  that  behind  a planar  wave.  Then  it  is 
readily  shown  that 


/P  A (r  2 - L2) 

C VBW  +1d‘L> 


(52) 


For  r » h it  follows  from  Equation  (52)  that 

0 


(P 

B~  C 


ftH- 


while  for  a cylindrical  confined  detonation 
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(53) 


fixed  distance  from  blast  center. 

As  indicated  by  Brode,et  al/11^  the  self-similar  blast  wave 
solution  used  here  provides  only  a crude  approximation  to  the  flow, 
and  with  recent  advances  in  numerical  computations,  more  accurate 
theoretical  solutions  have  become  available.  However,  as  was  shown 
by  the  illustrative  example,  r+  will  generally  be  small  compared  to 
the  fuel-cloud  radius  so  that  the  initiating  blast  will  not  have  a large 
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effect  on  the  total  impulse  which  is  generated.  The  propagation 

velocities  computed  using  the  zeroth  order  FAE  model  agreed 

remarkably  well  with  experiments  designed  to  simulate  a cylindrical 

FAE  with  infinite  confinement^. 

The  influence  of  side  relief  was  investigated  for  planar  waves 

but,  as  indicated  in  Part  H,  the  results  can  be  extended  to  cylindrical 

FAE's  for  sufficiently  large  values  of  the  wave  radius  r . The  in- 

s 

teraction  between  the  products  of  combustion  and  the  inert  boundary 
also  has  a self-similar  character  so  that  the  influence  of  side  relief 
and  cloud  height  can  also  be  represented  in  a remarkably  simple 
form. 

The  results  presented  here  make  it  possible  to  determine  the 
effects  of  fuel  properties,  cloud  geometry,  and  side  relief  upon 
impulse  with  a minimum  of  computational  effort. 


SECTION  III 


EXPERIMENTAL  RESEARCH 

A.  INTRODUCTION 

(8  9 

The  research  described  herein  is  part  of  an  ongoing  study^’  **  of 
cylindrical  shock  and  homogeneous  and  heterogeneous  detonation  waves, 
primarily  aimed  at  gaining  a better  understanding  of  the  fundamental 
aspects  involved  in  the  unconfined  explosion  of  fuel-oxidizer  mixtures. 
Some  of  the  early  results  of  these  studies  have  been  presented  at  meet- 
ings'*^; this  work  represents  extensions  made  since  the  previous 
annual  report.  Experimental  research  into  the  FAE  problem  has  been 
subdivided  into  two  major  phases;  studies  relating  to  the  initiation  of  self 
sustaining  homogeneous  and  heterogeneous  detonations;  and  studies  relat- 
ing to  the  breakup  and  ignition  of  liquid  fuel  drops.  The  former  phase  has 
naturally  received  most  of  the  emphasis.  The  latter  phase  was  under- 
taken to  examine  the  effects  upon  drop  time  history  of  the  trailing  rare- 
faction associated  with  the  blast  wave  and  the  attendent  reduction  in  gas 
pressure  and  temperature.  This  information  in  turn  will  assist  in  a 
more  complete  analytical  description  of  the  FAE. 

Regarding  the  first  phase,  experience  and  experimental  evidence 
substantiates  the  existence  of  four  regimes  of  explosion  which  may  occur 
upon  ignition  of  an  unconfined  fuel-air  cloud,  dependent  upon  the  intensity 
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of  the  initiating  source.  For  a low  energy  ignition  source,  such  as  a hot 
surface  or  small  open  flame,  the  cloud  will  likely  be  consumed  by  a normal 
deflagration  process  (although,  on  occasion,  they  have  been  observed  to 
accelerate,  presumably  due  to  turbulent  mixing  and  buoyancy  effects,  to 
full  fledged  detonations).  For  a somewhat  higher  and  rapid  energy  release 
(but  yet  sub -critical),  a blast  wave  followed  by  combustion  may  be  gene- 
rated which  then  rapidly  decays  to  a normal  deflagration  process.  For 

.*  * 

higher  rapid  energy  releases  (just-critical),  the  blast  wave  may  decay 
followed  by  a transition  stage,  at  some  distance  from  the  ignition  source, 
which  subsequently  leads  to  an  asymptotic  strengthening 
to  a CJ detonation.  For  yet  higher  rapid  energy  releases  (super -critical), 
the  blast  wave  decays  until  the  energy  contained  therein  is  of  the  same  order 
of  magnitude  as  that  due  to  combustion,  at  which  point  the  wave  begins  its 

. •.  ' k 

asymptotic  weakening  to  a CJ  detonation.  The  interest  of 
the  work  to  be  described  is  restricted  to  the  latter  three  cases  or  to  the 
case  of  a blast  wave  impacting  a single  stream  of  fuel  drops  in  an  oxidiz- 
ing atmosphere. 

The  experimental  facility  described  in  previous  report J®’ ^ was  used 
with  no  major  alterations  during  the  time  frame  reported  upon  herein. 
Following  the  acquisition  of  data  from  a series  of  blast  wave  shadowgraph 
photographs  , a slightly  improved  breech  was  developed,  and  more  efficient 
electric  ignitors  used.  An  extension  to  the  current  chamber  was  designed, 
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which  will  increase  the  chamber  length  to  nearly  6 feet,  in  order  to  accom- 
modate greater  run  times.  The  general  facility  capabilities  allow  initial 
drop  diameters  of  200  fi m < DQ  < 1800  fim,  with  any  non-corrosive 

oxidizing  atmosphere  of  0 < P < 2 atm  and  with  incident  blast  wave  Mach 

- ••  > • 

numbers  much  greater  than  common  Mach  numbers  of  detonation.  The 
Implementation  of  the  optical  system  used  with  the  chamber  is  discussed 
in  detail  in  Part  C,  where  it  received  the  greatest  use.  A final  word  is 

i , 

in  order  concerning  experimental  strong  blast  wave  data,  which  consumed 

• r ( 

much  time  to  compile  and  analyze  in  the  early  stages  of  this  work.  Typ- 
ical radius -time  data  obtained  in  the  sectored  chamber  can  be  non -dimen - 
slonallzed  and  displayed  as  shown  in  Fig.  16.  This  figure  demonstrates 
the  ability  of  the  chamber  facility  to  effectively  generate  blast  waves  of  a 

[highly  cylindrical  character;  something  which  has  been  stated  nearly  from 
the  start.  Figure  16  also  demonstrates  that  the  cylindrical  wave  charac- 
ter is  modeled  more  accurately  as  the  blast  source  energy  is  increased. 

B.  BLAST  WAVE  INITIATION  AND  PROPAGATION  OF  CYLINDRICAL 
DETONATIONS  ' * 

(1  4 9 16) 

Earlier  publications  ’ * ’ have  discussed  attempts  to  formulate 
satisfactory  models  for  the  blast  initiation  of  detonations  in  clouds  of  fuel 
in  gaseous  or  droplet  form,  and  pointed  out  their  Inherent  difficulties. 
Also  suggested  therein  was  a simplified  model  based  upon  a composite 
Of  similarity  solutions  for  the  strong  blast  and  CJ  detonation 


w 


(1 


(17) 

waves.  Bach,et  al v , in  an  attempt  to  overcome  some  of  the  difficulties 
alluded  to  above  and  still  possess  a model  of  sufficient  detail,  introduced 
a phenomenological  theory  of  initiation.  The  experimental  data  presented 
here  is  analyzed  and  discussed  in  light  of  the  advantages  offered  by  both 
of  these  approaches;  (1)  the  simplicity  of  the  composite  similarity  solution 
model,  and  ease  with  which  an  understanding  of  the  overall  process  is 
attained,  and  (2)  the  detail  of  the  phenomenological  model  with  which  an 
understanding  of  wave  transition  is  possible. 

The  essence  of  the  simplified  theory  is  used  as  a basis  for  a mathe- 
matical regression  model  developed  to  assist  in  the  analysis  of  the 
experimental  data.  First,  a brief  description  of  the  theory.  The  flow  is 
at  first  dominated  by  the  strong  blast  wave  with  transition  from  blast  to 
detonation  behavior  occurring  near  a critical  radius,  r+,  where  the  blast 
energy  and  energy  of  combustion  contained  within  r*  are  equal.  The 
complex  flow  in  this  transition  region  cannot  be  determined  analytically. 

In  the  simplified  theory  the  details  of  the  transition  region  are  ignored 
but  the  flow  is  represented  by  the  self-similar  solution  for  a strong  blast 
wave  for  r < r+  and  by  the  self-similar  detonation  solution  for  r > r+. 
Accordingly,  the  mathematical  regression  model  must  require  that: 

(1)  for  r < r*  the  wave  time  history  be  of  a second  order  dependency  on 
radial  displacement,  as  for  a strong  blast  wave,  and  (2)  for  r > r*  the 
wave  time  history  be  linearly  related  to  radial  displacement,  as  for  a 
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CJ  detonation.  Specifically,  the  model  makes  use  of  the 


method  of  LaGrange  multipliers  to  minimize  the  standard  error  of  a least 


squares  curve  fit  of  rough  experimental  time-radius  data  to  the  following: 


ti=(*0+*lri+<Vi  +el 


ri<  r* 


t,  = r,  + c 


1 ‘i  TCi 


ri>  r* 


The  problem  is  then  to  compute  the  constant  coefficients  a«,  a a 8 

O’  1’  V H0’ 

0J  and  the  critical  radius  r*,  in  order  to  minimize  the  standard  error 
subject  to  the  constraints  that  at  r[  = r+ 


a0+flf1r*+«2r* 


®1  +2a2r*  = ?! 


Confidence  limits  (standard  errors)  were  successfully  established  on  the 
evaluated  critical  radius  and  slope  of  the  linear  portion  of  the  model 
(i.  e,  detonation  velocity).  Typically,  detonation  velocity  is  evaluated  to  within 
2 percent  error,  but  the  error  on  critical  radius  often  exceeds  10  percent. 
Consequently,  in  practice,  the  detonation  velocity  derived  from  the  model 
is  used  unquestioned  but  the  critical  radius  value  is  tempered  by  manually 
tabulated  values  whose  error  is  of  the  order  of  the  model’s  value. 


Application  of  this  model  to  representative  experimental  data  will  be  discussed 
below.  Before  doing  this,  the  connection  between  the  simplified  theory,  Bach, 
et  al 's  phenomenological  theory,  and  the  experimental  data  must  first  be 
made  clear.  The  solid  lines  (1-8)  of  Figure  17  represent  the  essence  of  the 
latter  theory’s  behavior  on  an  Mg  versus  rg/ro  plot  where  rg  is  the  radius 
of  the  wave  and  rQ  is  the  explosion  length.  The  curves  essentially  repre- 
sent increasing  blast  source  energy.  The  model  recovers  the  two  limiting 
cases  of  a supercritical  energy  regime  (i.  e. , a monotonic  decay  to  a CJ 
wave,  curve  1)  for  very  large  source  energy  and  the  subcritical  energy 
regime  (i.  e. , the  blast  asymptotically  weakens  to  an  acoustic  wave,  curves 
5-8)  for  very  low  source  energy.  In  between  these  two  limits  (critical  energy 
regime)  the  blast  initially  decays  to  a sub-C  J velocity  and  then  asymptotically 
accelerates  to  its  eventual  CJ  condition  (curves  2-4) . The  model  predicts 
that  the  minimum  velocities  attained  in  the  critical  regime  depend  on  the 
magnitude  of  the  source  energy.  This  solution  suggests  that  unless  the 
initiation  energy  is  extremely  large,  a decay  below  the  C J condition  always 
occurs  and  the  eventual  approach  is  rather  slow.  By  comparison, the  simpli- 
fied theory  (curve  a)  predicts  an  immediate  initiation  of  the  C J condition 
upon  reaching  the  point  r*/rQ,  where 
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Figure  17.  Comparison  of  Blast  Initiation  Theories  , Simplified  and  Phenomenological 


r*  = (Eq/jt  pj  Q)l/2 

ro  = (Eo/2^yipi)1/2  (57) 

rJr0  = (2jr  Vj  Pj/Pi  Q)1^2 

This  ratio  is  a constant  for  a given  fuel  mixture  ratio  and  independent 

of  blast  source  energy.  Thus  on  an  M - r Jr  representation  all  blast 

s so 

initiated  detonations  collapse  to  a single  point  at  the  critical  radius  for 
the  simplified  theory.  It  is  now  clear  that  by  applying  the  regression  model 
to  experimental  data  which  presumably  follows  the  trends  of  say  curves  2 and 
3 in  Figure  17,  the  interpretation  suggested  by  curves  b and  c results.  Phys- 
ically it  is  apparent  that  a very  small  velocity  change  in  the  region 
0.  50  < rg/rQ  <0.  75  is  being  measured  and  consequently  being  referred 
to  as  a quasi-steady  sub-CJ  state.  Our  experiments  are 
limited  to  this  same  range  and  hence  a dependency  of  detonation  Mach 
number  upon  blast  source  energy  is  found. 

Figure  18  displays  in  transparent  fashion  typical  wave  behavior  of 
MAPP-air  detonations  for  9.  7 and  4.  3 percent  MAPP  by  volume,  respectively, 
as  a function  of  chamber  radius  and  source  energy.  Each  curve  represents 
a single  experimental  run  which  has  been  reduced  using  the  previously  dis- 
cussed regression  model.  The  experimental  detonation  Mach  numbers 
are  clearly  seen  to  be  monotonically  Increasing  with  initiation  energy. 

Runs  for  which  the  source  energy  was  sub-critical  are  represented  by  a 
second  order  least  squares  fit. 
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Figure  18a.  Experimental  Mach  No.  —Radius  Behavior,  9. 7 Percent  MAPP-Alr 
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Figure  18b.  Experimental  Mach  No.  —Radius  Behavior,  4. 3 Percent  MAPP-Alr 


General  features  displayed  by  Figure  18  are:  (1)  identifiable  threshold 
energies,  below  which  a quasi-steady  CJ  state  is  not  established;  (2) 
detonation  transition  distances  are  readily  shown,  beyond  which  the  calcu- 
lated detonation  Mach  number  is  achieved;  and  (3)  for  low  enough 
source  energy,  the  wave  process  proceeds  to  follow  the  pattern  of  a decay- 
ing reacting  blast  wave.  Additionally,  it  is  seen  that  curves  7 of  Figure  18a 
and  18b  both  suggest  the  establishment  of  detonation  wave  behavior,  yet  at  a 
clearly  sub-threshold  energy  level.  This  is  presumably  substantiated  by 
the  similar  occurrence  of  a stable  sub-CJ  wave  predicted 
by  the  phenomenological  theory  (Figure  17, curves  5 and 6) . A fundamental 
difference  between  Figures  18a  and  18b  is  suggested  by  the  fact  that  in  the 
former,  experimental  Mach  numbers,  Mg  , are  less  than  that  of  the  corres- 
ponding CJ  state,  (MCJ),  while  in  the  latter,  the  MCJ  is  straddled  by  the 
experimental  Mach  numbers.  The  reason  for  this  is  revealed  by  an  examin- 
ation of  the  orders  of  magnitude  of  blast  source  energy,  E,  used  in  experi- 
ments at  a given  mixture  ratio,  relative  to  the  corresponding  critical  energy, 

E . In  the  9.7  percent  MAPP  case  the  ratio  E/E^  never  exceeded  1.2,  while  in  the 

4.  3percent  MAPP  case  this  ratio  approached  9.  5.  The  situation  whereby 

M > M„  , typifies  the  supercritical  energy  regime  for  wave  propagation 
exp  v/j 

and  the  above  results  suggest  that  this  region  is  established  as  E/E  ^ — 10.0. 

This  statement  is  more  easily  demonstrated  by  examination  of  the  actual 
experimental  data  given  in  Table  II  and  displayed  graphically  in  Figure  19. 


Extrapolation  of  existing  data. 


gation  was  experimentally  observed  is  easily  seen.  It  is  of  interest  to  examine 


this  line  of  thought  in  light  of  MAPP-air,  as  well  as  available  two-phase 
detonation  Mach  number  data.  Plotted  in  Figure  20  is  detonation  wave  Mach 
number  against  volumetric  fuel  concentration  for  curves  of  constant  blast 
source  energy.  The  solid  lines  represent  experimental  MAPP-air  data, 
with  available  kerosene-air  two-phase  data  superimposed.  Also  given  is  a 
curve  of  MCJ  for  MAPP-air  with  the  corresponding  data  given  for  kerosene- 


I 


air.  A line  of  may  be  constructed  by  the  intersection  of  the  appropriate 

curve  of  constant  E with  the  appropriate  fuel  concentration.  Hence,  the  line 

0f  Ecrit  f0r  MAPP-air  data  is  shown,  again  with  the  corresponding  kerosene- 

air  point  also  shown.  The  importance  of  this  figure  is  revealed  when  it  is 

observed  that  for  the  MAPP  concentrations  which  experienced  M > M 

exp  — CJ 

their  non-dimensional  energies  were  2. 8 < E/Ecrit  < 9.  3.  Correspondingly, 

at  the  given  kerosene  concentration  was  nearly  equivalent  to  M with 

E/Ecrit  “ 2'5  ' tt  is  thus  apparent  that  for  E/EcrU  approximately  greater 

than  3.0,  wave  propagation  characteristics  associated  with  the  super-critical 

energy  regime  can  be  established  and  studied.  The  interval  1.0  < E/E  < q n 

' crit  * u 

then  approximately  specifies  the  critical  energy  regime,  which  does  indeed 
seem  to  be  very  narrow.  This  discussion  appears  to  break  down  as  the  fuel 
detonabillty  limits  are  approached,  as  shown  in  Figure  20.  In  lean  mixtures 
nearly  all  non-dimensional  energies  suggest  supercritical  wave  behavior. 
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Comparison  of  Experimental  Mach  Number  Behavior  with  Theory  Plotted 
Against  Fuel  Concentration  with  Varying  Blast  Source  Energy. 
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This  Is  presumably  due  to  incomplete  chemical  reactions  as  the  detonation 


limits  are  approached,  which  is  easily  responsible  for  M deviating 

CaU 


from  The  overall  detonation  wave  speed  behavior  observed  in  MAPP- 

LJ 


air  experiments  is  seen  in  Figure  21.  Plotted  is  the  range  of  experimental 


Mach  numbers,  M , obtained.  As  was  already  pointed  out  in  the  previous 

exp 


discussions,  this  range  is  indeed  not  to  be  interpreted  as  a range  for  experi- 
mental error,  but  rather  a function  range  within  which 

M = M(E)  (58) 


It  is  foreseen  that  through  the  use  of  spark  schlieren  data,  wave  structure 
may  be  analyzed  as  source  energy  is  varied  throughout  the  multiple  energy 
regimes.  Perhaps  then  a more  precise  estimate  of  non-dimensional  energy 
intervals  may  be  established. 

Experimental  transition  distances  when  non-dimensionalized  by  explo- 
sion length,  where 


r*/r0  = 2/Md  (y22  - 1)1/2 


(59) 


for  cylindrical  wave  behavior,  compare  favorably  with  those  computed 
based  upon  theoretical  detonation  Mach  number.  Figure  22  is  a plot  of 
non-dimensional  critical  radius  against  fuel  concentration  with  experimental 
MAPP-air  data  superimposed.  Also  shewn  is  the  available  kerosene-air 
data  for  comparison,  with  attendant  uncertainty  displayed  in  all  data. 
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Some  deviation  between  experiment  and  theory  is  expected  in  view  of  the 
observation  that  exact  CJ  detonation  is  usually  not  attained.  Further, 
even  the  idealized  theory  of  strong  blast  wave  behavior  up  to  r*  and  CJ 
detonation  beyond  r*  does  not  predict  the  attainment  of  the  CJ  Mach  number 
exactly  at  r*. 

Wave  pressure  histories  at  discrete  radial  locations  have  also  been 
determined  experimentally.  Shown  in  Figure  23  are  typical  pressure 
histories  at  r = 19.  5 inches  and  r = 30.0  inches  for  strong  blast  waves 
and  gaseous  and  two-phase  detonations,  at  two  source  energies.  Similar 
volumetric  fuel  concentrations  were  chosen  for  this  comparison  with  4.  3 
percent  MAPP  for  the  gaseous  detonation  data,  and  4. 1 percent  kerosene  for 
the  two-phase  detonation  data.  The  higher  source  energy  initiated  a detona- 
tion in  both  of  these  cases,  while  the  lower  energy  did  not.  Examination  of 
such  data  clearly  displays  several  fundamental  features:  (1)  identifiable 
difference  between  pressure  levels  of  detonating  and  non -detonating  cases, 
(due  to  wave  Instabilities  this  difference  is  not  as  distinctive  at  E .,); 

(2)  higher  sustained  pressure  for  gas  phase  over  two-phase  detonation; 

(3)  approximate  location  of  transition  distance  is  identifiable  in  the  traces; 
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Figure  23.  Typical  Wave  Pressure  Histories  at  Radial  Distances  of  R = 19.5 
and  R = 30.0  Inches  for  Blast  Source  Energies  of  2.0  Grams  and  0.0  Grams 
Detasheet.  (Upper  Traces  are  for  R = 19.5,  50  ^xsec/cm,  Lower  Traces 
are  for  R = 30.0  Inches,  20  fisec/cm.) 
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whereas  the  two-phase  mixture  detonated  between  the  pressure  transducers 
(downstream  pressure  higher  than  upstream),  the  gas  phase  mixture  detonated 
prior  to  the  first  pressure  transducer;  and  (4)  clear  differences  in 
pressure  from  the  pure  blast  wave  case,  even  for  the  non -detonating  cases. 

C.  BLAST  WAVE  BREAKUP  AND  IGNITION  OF  LIQUID  FUEL  DROPS 
The  final  phase  of  the  current  experimental  exploration  into  the  funda- 
mentals of  blast  initiation  of  an  unconfined  fuel-cloud  involves  the  study  of 
individual  drop  dynamics.  The  test  conditions  examined  during  this  phase 
included  two  fuels,  normal -propyl  nitrate  and  decane  with  drop  sizes  of 
384fi  and  768ju,  for  three  ambient  oxidizers  of  100%  00,  50%  Oo-50%  Nnand 
air  at  1 atmosphere,  and  for  incident  blast  wave  strengths  1.  5 < M < 6. 

The  essential  findings  are  enumerated  in  detail  in  Table  III  in  the  form  of 
essential  test  conditions,  drop  breakup,  stripping,  and  ignition  data,  as  well 
as  overall  drop  dynamic  conditions  initially  and  at  breakup.  Such  data  was 
obtained  from  a computer  program  developed  to  analyze  in  detail,  digitized 
representations  of  streak  schlieren  records  obtained  during  test  runs. 

The  optical  system  used  Ls  shown  in  Figure  24.  Shown  in  Figure  25  are  typical 
streak  schlieren  records  of  the  breakup  and  ignition  of  768jn  drops  of  n- 
propyl  nitrate  in  air  and  100%  C>2  at  standard  conditions,  for  two  different 
source  energies.  In  the  records  the  time  axis  is  running  horizontally,  for 
a right  running  shock,  with  the  distance  axis  running  from  bottom  to  top. 

The  three  horizontal  lines  are  reference  wires  at  1 inch  intervals,  with  which 


TABLE  in.  EXPERIMENTAL  RESULTS  OF  INDIVI DUAL  DROP  DYNAMICS  STUDY. 


* 


i I*  ill  ii’»  a \ 5 *«  .*$§  j | 

HHI9II  l n 11s  m 13  i 5 II  iU  I l 

ssl=E»s2  i s*  jjjji  if*  is  i I i!§  iss  i i 

MMMM  5 58  §11  58s  58  5 i it!  Hi  1 i 


gsg  z s 

808  S 2 


* ! 

I EsE.SSS. 

0 

-5 

M — 

• 8 

o 

5 S**|*=c4 

i 

s« 

U- 

sii 

2 

8 

ss* 

aii 

f 

iisiiilg 

c 

X 

li 

555 

ill 

if 

5 

8 

iii 

III 

i '* 

55851115  i 

55 

ill 

555  15 

8 

S 

§ 

ill 

115 

•“  'o 

K 

aismgi 

1 

IS  iii 

ill 

15  ! 

i 

511 

sss 

288 

888. .So* 

_ 

o O 

* 

► i 

S *a=*2*3= 

8 

a* 

EE! 

iii 

s* 

5 

3 

sis 

888 

.»  i 

i ::!:«!! 

e 

: 

?E: 

^ : 8 
o'  o 

:*  = 
o o 

! :« 

8 

: 

553 

: as 

> $ 

p . . 

*F. 

. , s 

. 

S 

sggc 

S3 

* * Js 

sfe- 

8 ' 3 

* is 

' ' * 

s 

* 

ssst; 

j* 

55...EL 

o 

0 

.j 

§ 

0 

g 

0 

o 

o' 

boo 

- o 

o 

• 

b 

b 

bob 

.*  ; 

e 

o 

o o 8 

i s 

o 

o 

1 

o 

° 

boo 

o o 

° 

b 

° 

° 

ooo 

SIIsi!*S 

o 

8 

ASS 

1 2 

8 

8 

? 

i 

S;  s 

•bbbbbbb 

• 

bob 

“ • 

• 

b 

< j 

Slililis 

E 

s 

sag 

3 3 

8 

8 

i 

8 

32s 

-*  ! 

° 

b 

ooo 

o b 

° 

° 

b 

b 

ooo 

. 1 

! i|!?*SS2 

E 

5 

ESI 

= 8 

X 

a 

M 

S 

233 

- 

• 

bob 

8 

b 

b 

b 

b 

bbb 

*E!ESI*i 

i 

Si 

ale 

as* 

ES 

3 

s 

5 8 

23: 

2 b •’  b b b b b 

4 

bb 

4*4 

2»‘» 

8 

82 

328 

, 88288288 

s 

as 

3 S 8 

*88 

PS 

8 

8 

SS 

3 8S 

! bbbbbbbb 

• 

b o' 

bob 

bob 

o^ 

b 

b b 

ElliiSEi 

8 

If 

ih 

EEl 

§3 

E 

l 

33 

§S8 

• 

bb 

b 4b 

bbb 

wb 

b 

bb 

<1 

■itrtiii 

8 

88 

5 3 8 

688 

" W 

s 

8 

xs 

X 8 8 

444 

bbb 

b*4 

" 

•* 

b b 

«.o 

1 

1 88888888 

8 

88 

888 

12? 

ss 

a 

«• 

33 

838 

j 

! 4s444c*4 

4 

44 

*44 

838 

is 

4 

b 

ss 

S8X 

4 

S' 

Islf IEIII* 

EiSitf 

HEHiHS 

iSHS 

*»Si 

/ 

444444 

bbbbbbbbb 

•i  ti  ti  ti  m 

b 4 b < 

/■ 

12*21221 

44.44444 

1 

b 

El 

44 

HE 

bbb 

iSS 

bbb 

as 

b b 

i 

b 

i 

!l! 

bbb 

21* 

bbb 

I I 

8 8 8888* 

■88888 

* ***E*tt 

88888 

888888  8 . 

b b 4««4b 

444444 

b b b b b b b b 

b b b b b 

O-bbbb  • 4 

J 

1 

1 

1 

<nr 

» 

li 

* *“  3 

I 

J 

!>*  ]i 


ssssic  tssxsssas  ?:»:*  «**sss  s 

73 


TABLE  m.  EXPERIMENTAL  RESULTS  OF  INDIVIDUAL  DROP  DYNAMICS  STUDY 

(Concluded). 
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Figure  25a.  Typical  Streak  Schlleren  Record  of  768^  Drops  of 
n-Propyl  Nitrate,  Incident  Mach  Number  of  2.0  in  Air. 
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Figure  25b.  Typical  Streak  Schlieren  Record  of  768/jl  Drops  of 
n-Propyl  Nitrate,  Incident  Mach  Number  of  3.  5 in  Og* 


76 


a nominal  time  of  50  nsec  along  the  horizontal  can  be  associated.  The 
fuel  drop  enters  the  field  of  view  from  the  left,  it  is  impacted  by  the  shock 
entering  from  the  bottom,  the  drop  is  subsequently  accelerated  in  the 
downstream  direction,  simultaneously  shattered  and,  provided  proper 
conditions  exist,  ignited.  These  photographs  clearly  display  the  changing 
local  dynamic  conditions  the  drop  experiences  during  its  breakup  and  igni- 
tion history.  The  source  energy  of  Figure  25b  was  sufficient  to  cause 
ignition.  In  this  case  the  observed  ignition  time  is  approximately  45  jasec 
for  an  incident  Mach  number  of  3.  5.  This  value  is  substantially  higher 
than  similar  data  reported  by  Lu  and  Slagg^8}  Presumably,  this  difference 
is  due  to  the  trailing  rarefaction  associated  with  the  blast  wave  and  the 
attendant  reduction  in  gas  pressure  and  temperature.  This  work  is  unique 
with  regard  to  this  key  point,  since  all  previous  drop  breakup  and  ignition 
studies  were  confined  to  the  simple  case  of  an  incident  plane  shock  wave, 
with  which,  during  the  testing  interval,  no  such  reductions  are  associated. 

A measure  of  the  magnitude  of  reduction  in  local  drop  dynamic  conditions 
is  given  by  the  ratios  q^/q^  and  He^/Re^,  which  compare  the  dynamic  pres- 
sure and  Reynold's  number  lor  the  drop  initial-  and  breakup -conditions. 

The  data  tends  to  suggest  that  in  general  this  reduction  is  quite  severe, 
with  typically 


A final  reference  to  Figure  25b  will  reveal  the  presence  of  combustion 
products  carried  by  the  convective  gas  flow  a distance  of  approximately 
65  iis  (at  the  drop  location)  behind  the  incident  shock.  These  products 
originate  at  the  blast  source  and  serve  to  point  out  the  important  difference 
between  explosive/ignitor  and  spark  discharge  or  even  exploding  wire 
blast  wave  energy  sources.  This  difference  of  course  being  the  addition  of 
mass  to  the  flowing  system  upon  generation  of  the  blast  wave.  Practical 
experimental  considerations  led  to  the  observation  that  for  too  high  a blast 
source  energy  (approximately  2. 5 grams  of  Detashee$,  all  details  of  the 
combustion  process  were  literally  obscured  from  view.  Some  courses 
of  action  remain  open  for  possible  alleviation  of  this  problem  if  further 
tests  are  warranted.  Further  detailed  analyses  of  the  data  reported  herein 

-m 

are  planned. 
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